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^ • Abstract 



Moments of the renormalized i?-meson shape function provide a natural way to define 
' short-distance, running heavy-quark parameters such as the 6-quark mass and kinetic 

^ energy. These parameters are particularly well suited for studies of inclusive decay 

distributions. The definitions of nib and //^ in this "shape- function scheme" are derived 
to two-loop order. Using previous determinations of heavy-quark parameters in other 
schemes, we find m6(///) = (4.63 ± 0.08) GeV and nUi^f ) = (0.15 ± 0.07) GeV^ at a 
reference scale fJ-j = 1.5 GeV. 



Introduction. The past decade has seen a revolution in the precision with which funda- 
mental physics can be probed using measurements of the decay properties of b quarks. For 
instance, the element \Vcb\ of the quark mixing matrix is now known with a precision of 2%, 
which is close to the accuracy of our knowledge of the Cabibbo angle PJ. Measurements of 
rare decay processes such as i? ^ Xsj impose stringent constraints on model building. The 
6-quark mass has been determined with an uncertainty of about 60MeV P], which is much 
less than the QCD scale. To achieve such precision requires that heavy-quark parameters 
can be defined unambiguously using short- distance techniques. In particular, these definitions 
should not refer to the notion of on-shell quark states, which would introduce uncontrollable 
uncertainties ("renormalon ambiguities") that far exceed the experimental precision achieved 
at the B factories. 

The question of a short-distance definition of the 6-quark mass has received much attention. 
While an ad hoc subtraction scheme such as MS in principle defines the quark mass in an 
unambiguous way, such a definition is not appropriate for the discussion of B decays, where 
the typical scales are often significantly below mi,. A more fruitful concept is that of a low-scale 
subtracted quark mass |2], which is based on the idea that non-perturbative contributions to 
the heavy-quark pole mass can be subtracted by making contact to some physical observable. 
The result is an expression m;,(yUj), which differs from the pole mass by an amount proportional 
to a subtraction scale /i/ = few x Aqcd- 

Several examples of such low-scale subtracted quark masses have been discussed in the 
literature. In the potential-subtraction scheme, long-distance contributions to the pole mass 
are subtracted by relating it to the static potential between two heavy quarks jS] . Similarly, in 
the T(15') scheme the quark mass is related to the mass of the lowest-lying bottomonium res- 
onance Both schemes are well suited to study heavy-quark systems in the non-relativistic 
regime, such as bb spectroscopy or heavy-quark production near threshold. A quark-mass 
definition more closely related to the non-perturbative physics probed in i?-meson decays is 
provided by the kinetic scheme [3121111; in which the non-perturbative subtraction is accom- 
plished with the help of heavy-quark sum rules [71. Because these sum rules constrain the 
properties of i?-meson form factors in the "small- velocity limit", the kinetic scheme is well 
suited for heavy-quark expansions applied to B decays into charm particles. 

Inclusive B decays into final states consisting of only light hadrons, such as B Xs'y or 
B — > Xu 1 1^, probe yet different aspects of non-perturbative physics. The bound-state effects 
relevant in these processes are encoded in 5-meson shape functions defined in terms of the 
forward matrix elements of non-local string operators on the light cone |H1 El • A sensible 
definition of heavy-quark parameters for such processes should incorporate the bulk properties 
of the leading-order shape function S{uj). Schematically, inclusive decay spectra are given in 
terms of convolution integrals of the form 

I du f{mT'^ + u) S{u) ^ f{ml°'^ + {u)) ^ f{mf) , (1) 

where / is some distribution, and in the last step we have defined the "shape-function mass" by 
adding the average value of u to the pole mass. More precisely, the basis of the shape-function 
scheme proposed in is to obtain short-distance definitions of all heavy-quark parameters 
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that enter the description of inclusive decay rates via the moments 

/oo 
duu^ S{uj,fi) (2) 

of the renormahzed i?-meson shape function S{u!, fi), which has support over the range — oo < 
uj < rriB — iTLh- At tree level, the shape function is normalized to unity, its first moment 
vanishes (in the pole scheme), and higher moments are given in terms of local operator matrix 
elements in heavy-quark effective theory, e.g. Mg''*^'^ = — Ai/3 |H1E], where Ai is the kinetic- 
energy parameter defined in the pole scheme ^21- While it is well known that the 5-quark pole 
mass suffers from a renormalon ambiguity, the same is true for other heavy-quark parameters 
such as Ai [13], even though this may not be obvious from low-order perturbative calculations 
|14j . Beyond the tree approximation, the moments Mjv(yUj,/i) depend on the renormalization 
scale fj, and the lower cutoff fif applied to the integral in (j2|), and it is natural to use them 
to define a set of running heavy-quark parameters order by order in perturbation theory. In 
fact, for /i/ ^ Aqcd these moments can be calculated using an operator product expansion, in 
which the scale Hf plays the role of a hard Wilsonian cutoff. The moment integrals can then 
be matched onto a series of forward S-meson matrix elements of local operators. The Wilson 
coefficients in this matching can be obtained using perturbation theory with free quark and 
gluon states. In particular, if the matching is performed in the parton model with on-shell 
heavy-quark states {v ■ k = 0, where k = pf, — mh v is the residual momentum, and v denotes 
the 5-meson velocity), then the local operator matrix elements are given by their tree-level 
expressions, and the matching becomes trivial once the moments (j21) are computed to a given 
order in as- 



Two-loop relations in the shape-function scheme. The parton-model expression for 
the renormahzed shape function was derived at one-loop order in pT^ llSj. The result involves 
so-called star distributions JH]; which are generalizations of plus distributions appropriate for 
test functions defined on an unbounded interval. This complication can be circumvented by 
introducing the function 



s(L,/i) 



duj S, 



parton 



In 



fif + n ■ k 



(3) 



whose dependence on the cutoff fif enters only through the logarithm L, while /i dependence 
also resides in the coupling as{fi)- Here n is a light-like vector (n^ = 0, v - n = 1), which enters 
in the definition of the shape function as the matrix element of a non-local light-cone string 
operator |9|. The shape-function moments can now be obtained as 



M^(/i^,/i) = (-1) 



N 



In 



a + n ■ k 
/i 



(4) 



At the end of the calculation one must expand the result in powers of n-k and identify n-k ^ 0, 
{n-kf -Ai/3, etc. [II]. 

The function s{L,fi) obeys an integro-differential renormalization-group equation, which 
can be derived starting from the evolution equation for the shape function obtained in [TT] ll7j. 
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We find 



s{L, fi) = 2 rcusp(«s) L - 7(as) s{L, /i) 



+ 2T 



1 dz 



z 



s{L + ln(l — z), jj) — s{L, /i) 



(5) 



where as = o:s{^), and Fcusp, 7 are anomalous dimensions. This equation can be solved order 
by order in perturbation theory with an ansatz of the form 



n=l 



Ait 



(6) 



The evolution equation (jSJ allows us to express all coefficients of logarithms in terms of the 
perturbative expansion coefficients of the anomalous dimensions and f3 function, defined as 



n=0 

dlnfi 



as 



n+l 



n=0 



n+l 



n=0 



a. 



n+l 



(7) 



At two-loop order, we obtain 



= 1 + 



+ 



Air 

. Att 



c?^+27oL-roL2 



) + (24^^(70 - /?o) + 271 + ^ ro7o + 4C3 J L 



+ 1^270(70 - Po) - ci'^Fo -T^--T■',jL'+[^^(3o- 2jo) ^oL' + - L 

The relevant expansion coefficients of the f3 function and cusp anomalous dimension r^usp 
are (in the MS renormalization scheme) 



^2 t4 



R 2 



T^ = ACf^ 



'268 _47r2^ 
~9 3~ 



CA--nf 



(9) 



The two-loop coefficient of the anomalous dimension 7 has been calculated in JH]- We have 
found some mistakes in the translation of the results for the two-loop graphs into the expression 
for the anomalous dimension. The corrected result is |17| 1201 



7o 



-2Cf 



11 =C, 



74 vr 



+ — - I8C3 + Ca + — + 



TT 



27 18 



27 9 



(10) 



where k = under the assumption that the two-loop diagrams themselves were evaluated 
correctly in pSJ. However, there is reason to believe that there might be an additional error 
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in that paper, giving rise to a non-zero value k = 4/3 [21], which we adopt in our numerical 
analysis. The expression for 71 should be checked with an independent calculation. Finally, 



the coefficient c. 



-7r^/6) Cp was computed in [Ullin], while the coefficient Cg is presently 
unknown. This coefficient does not enter the definitions of heavy-quark parameters at two-loop 
order in the shape-function scheme. 

Using the explicit two-loop expression (jSJ for the function s{L,fi), it is straightforward to 
compute expressions for the shape-function moments (jlj) in the pole scheme. Closed expres- 
sions valid at one-loop order are given in JTT . One finds that the first moment Mi{fif,fi) no 
longer vanishes beyond the tree approximation but receives a perturbative correction propor- 
tional to fijasifi). The key idea behind the shape-function scheme is to redefine the 6-quark 
mass, ml°^^ = mh{fLf,fi) + 6m, in such a way that after the redefinition the first moment 
vanishes. Technically, this is implemented by introducing a residual mass term 6m for the 
heavy quark |22j. All that changes in the expressions for the moments is that the residual 
momentum k gets replaced hj k + 6m v, and hence n ■ k ^ n ■ k + 6m. We then choose 
6m such that Mi(/ij-,/i) = 0, thereby defining a mass renormalization scheme order by order 
in perturbation theory. Next, we define the kinetic-energy parameter via the ratio 

3M2(/i/, /u)/Mo(/U/, yu) in the new scheme. This quantity differs from the pole-scheme param- 
eter — Ai by terms of order fijasifi). We then solve for m^°'^ and Ai in terms of the new 
parameters mbifJ^j, /i) and fil{fif,fi). To two-loop order, the results of these manipulations are 



IT 



+ 



3/i/ 



TT 



l-21n^ + ^ A:, 



21n^ + 



/X vr 



A;2(/i/,/i) 



-Ai = 



1 + f-i - 3 In + ^ k,{^^„ 

TT \ 2 jj, TT 



2 CFas{^i) 



TT 



's^n^ + ^k,i^^J,^^) 

H TT 



where the one-loop terms agree with ^T], while the two- loop corrections are encoded in the 
coefficient functions 
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TT 



27 
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24 



— TT 
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7r 
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4 / 



4Ci.ln2^ 
/i 



(12) 



Equations (jlip and (jl2j) are the main results of this work. They allow us to consistently 
implement the shape-function scheme with two-loop accuracy in perturbative calculations of 
heavy-quark processes. They also enable us to connect the shape-function scheme with any 
other short-distance scheme used to define heavy-quark parameters. 

For practical applications, it is often inconvenient to use heavy-quark parameters defined 
as a function of two scales, and we will therefore adopt the special choice /i = /x/ as our 
default, denoting mf,(/ij) = mb(/ij,/ij) and = /i/). Then the relations (fTTj) 

simplify, since all logarithms vanish. For = 3, we obtain 



pole 



mb{fif)+fif 
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TT 
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(13) 



Connections with other schemes. To the best of our knowledge, the only other scheme 
in which both and /i^ are defined in a consistent way (without renormalon ambiguities) is 
the kinetic scheme j2] . The relevant two- loop relations in that scheme are El 



pole 

ml 



4 CpOLsi^iif) 
mbi^^f) + ^ 



1 + 



IT 



13 
12 



TT 
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^ fx} CFas{2fif) r ^ as{2fif) 

2mb{fif) TT \ TT 

^ 2 CFas{2iif) ^ a.(2/i/) 
-Ai = fi^lfif) - fif jl + 



13^ /13 7r2\ 

12^°+ 1,12 -yj^^ 

13^ /13 7r2\ 
12^°+ 12~Yr^ 



(14) 



The fact that the 6-quark mass receives a l/ruh correction in the kinetic scheme leads to a 
mixing between different terms in the l/rrif, expansion. No such mixing occurs in the shape- 
function scheme. Other approaches such as the potential-subtraction or T(IS') schemes only 
redefine the 6-quark mass but ignore the problem of infrared ambiguities of other heavy-quark 
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parameters. The two-loop relation between the potential-subtracted quark mass and the pole 
mass reads |2] 



TT 



1 + 



a 



11 2 



(15) 



At one-loop order, the potential-subtracted quark mass coincides with the mass defined (at 
the same scale fif) in the shape-function scheme. The next-to-leading order relation between 
the pole mass and the 6-quark mass defined in the T(IS') scheme is 



pole 

mr - 



8 



■ m 



TC 



L + 



11 



)/3o- 



(16) 



where L = Inlfi/Cpasifi) 1^1°^^]. The peculiar counting of powers of as{^) in the T(15') scheme 
becomes more transparent if we rewrite this relation in the form 



1 + 



TT 



TT 



In 



^/i/(^) 6 / ° 3 ^ 



(17) 



where 



TT 



pole 



plays the role of the subtraction scale. Note that the T(IS') mass is scale independent, contrary 
to the other low-scale subtracted quark masses considered above. On the other hand, the 
subtraction point yU/(yu) is itself a function of the renormalization scale fi. Based on (fTB|) one 
would suspect that a proper choice of /i is such that L ^ 0, which yields yU ~ 2 GeV. However, 
relation (fTH|) suggests that the appropriate physical scale in the T(15') scheme may be lower, 
in accordance with arguments presented in 



Numerical results. Estimates of the 6-quark mass and kinetic energy in the shape-function 
scheme can be obtained using various sources of phenomenological information derived in other 
subtraction schemes. For this purpose we adopt the scale choice fif = fi = 1.5 GeV in the 
shape-function scheme. The running of the heavy-quark parameters will be studied later. As 
mentioned above, we use k = 4/3. 

The values of rrih and /i^ defined in the kinetic scheme at a reference scale /i/ = 1 GeV have 
recently be determined from a global fit to experimental data on moments of various B —>■ 

/ decay distributions P, yielding m;, = (4.611±0.068) GeV and /x^ = (0.447±0.053) GeV^ 
(with a mild anti-correlation of errors, c = —0.4), in good agreement with theoretical expec- 
tations Using these results, we compute the pole scheme parameters from ()14|) and then 
solve for the parameters in the shape-function scheme using (fT^ . We obtain 

mfe(1.5 GeV) = (4.63 ± 0.07) GeV , /i^ (1.5 GeV) = (0.15 ± 0.06) GeV^ , (19) 

where the errors refiect the uncertainties in the kinetic-scheme parameters only, and the cor- 
relation is the same as in the kinetic scheme. Note that we solve the relation for /i^(/i/) in 
(fT^ exactly, rather than inverting it to O(a^). (In the latter case the result for nib would 
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Figure 1: Dependence of the heavy-quark parameters in the shape-function scheme 
on the scale v at which the coupling constant is evaluated. The lines refer to the 
translation of results obtained in the kinetic scheme (solid), the potential-subtraction 
scheme (dashed), and the T(IS') scheme (dashed-dotted). Central values for the 
heavy-quark parameters in these schemes are used as input. The points show the 
values in for comparison. 



stay unchanged, whereas the value of yU^ would increase to (0.18 ± 0.06) GeV^, which is in- 
side the error range.) For reference, the corresponding results obtained at one-loop order are 
mb{l.h GeV) = (4.57 ± 0.07) GeV and /x^ (1.5 GeV) = (0.34 ± 0.06) GeV^. The significant im- 
pact of two-loop corrections in the case of /i^ is due to the terms proportional to /ij in the 
relation for Ai in (fT^ . which arise first at 0{a1). 

The potential-subtracted mass at the scale /i/ = 2 GeV has been determined from an 
analysis of the bb production cross section and the mass of the T(IS') resonance to be mj, = 
(4.59 ± 0.08) GeV |2lj. Using relations (fT3|l and (fT5|) . as well as the value of /x^ given in 
(fnH) . we may use this information to obtain mfe(1.5GeV) = (4.66 ± 0.08) GeV for the h- 
quark mass in the shape- function scheme, in excellent agreement with (jl9p . We may also 
use the value mb = (4.68 ± 0.03) GeV obtained in the T(IS') scheme from a fit to moments 
of B —>■ Xcliy and B Xsj decay distributions [2^1 , even though we do not trust the 
very small error assigned in this analysis. Taking /i = 1.5 GeV for the scale in |T6|) . we find 
mb(1.5 GeV) = (4.61 ± 0.03) GeV for the 6-quark mass in the shape-function scheme, which is 
again consistent with ()19j) . 

The numbers presented above are affected by additional theoretical uncertainties inherent 
in the translation from one scheme to another. A way of assessing these uncertainties is to 
evaluate the coupling constant ««(//) in the various relations at a common scale u different 
from the individual "natural" scales. This can be done using 



a^(/i) = a^(i/) 



l + ^^ln^ + ...' 

2 TT /i 



(20) 



While the results are formally independent of the choice of u, the variations due to truncation 
effects may serve as an indicator of the impact of higher-order effects. Varying between 1.5 
and 2 GeV so as to cover the range of "natural" scales in the various schemes yields the results 
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shown in Figure ^ The perturbative uncertainties indicated by this analysis are somewhat 
smaller than the parameter uncertainties shown in (jl9p . Combining the two sources of errors, 
we quote our final results as 

mb(1.5 GeV) = (4.63 ± 0.08) GeV , /i^(1.5 GeV) = (0.15 ± 0.07) GeV^ . (21) 

Ultimately, the shape-function parameters should be determined directly from experimental 
data, without a detour through another subtraction scheme. Moments of the photon energy 
spectrum in B ^ Xs'j decay are ideally suited for this purpose and will allow an extraction 
of rub with far better precision j27j . 

Given values of the heavy-quark parameters mb and f^l defined in the shape-function 
scheme at some reference scale /x/, it is straightforward to solve the relations (fT!?|) to find 
the corresponding running parameters at a different scale. There is nothing unphysical about 
this scale dependence, which indeed is an essential feature of our approach. For instance, 
we obtain mfe(l.lGeV) = (4.64 ± 0.08) GeV and /i^(l.l GeV) = (0.18 ± 0.07) GeV^ at a 
lower scale, which is appropriate for applications the B —>■ Xg'j photon spectrum ^Tj ITT] , 
and 7716(2 GeV) = (4.60 ± 0.08) GeV and /i2(2GeV) = (0.11 ± 0.07) GeV^ at a higher scale, 
which could be used to analyse moments oi B ^ Xcl i' decay distributions. More important 
evolution effects are encountered when the scales /i/ and /i are taken differently, in which 
case the relations (fTT|) must be used to compute the effects of scale changes. For example, 
with fxf = 1.1 GeV and /i = 1.5 GeV we get 777^(1.1 GeV, 1.5 GeV) = (4.58 ± 0.08) GeV and 
liKl.l GeV, 1.5 GeV) = (0.36 ± 0.06) GeV^. 

Conclusions. In summary, we have derived the two-loop relations for the 6-quark mass 
rrib and kinetic-energy parameter /i^ in the shape-function scheme introduced in ^T]. These 
relations are important for performing high-precision studies of -B-decay observables. We have 
determined values for mf,(/7/) and ^K^j) at a reference scale /i/ = 1.5 GeV by using as input 
the corresponding values of heavy-quark parameters in other subtraction schemes. In the 
future, the shape-function parameters should be determined directly from a fit to moments of 
the B — > Xg'j photon spectrum. 
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Note added. After this work was submitted for publication, an Erratum to jTHI appeared, 
in which the authors confirm the result (fTUj) with tc = 4/3. 
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